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Fulde-Ferrell-Larkin-Ovchinnikov Superconducting State in CeCoIn 5 
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We report specific heat measurements of the heavy fermion superconductor CeCoIn 5 in the vicinity 
of the superconducting critical field H C 2, with magnetic field in the [110], [100], and [001] directions, 
and at temperatures down to 50 mK. The superconducting phase transition changes from second to 
first order for field above 10 T for H || [110] and H || [100]. In the same range of magnetic field we 
observe a second specific heat anomaly within the superconducting state. We interpret this anomaly 
as a signature of a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) inhomogeneous superconducting state. 
We obtain similar results for H | [001], with FFLO state occupying a smaller part of the phase 
diagram. 

PACS numbers: 74.70.Tx, 71.27.+a, 74.25.Fy, 75.40.Cx 
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In the early 1960's, following the success of the BCS 
theory of superconductivity, Fulde and Ferrell pj and 
Larkin and Ovchinnikov [2| developed theories of inho- 
mogeneous superconducting states. At the core of FFLO 
theory lie competing interactions of a very basic nature. 
One is the interaction of the spin of the electron with 
magnetic field and the other is the energy of the super- 
conducting coupling of electrons into Cooper pairs, or 
the condensation energy. In the normal state the elec- 
trons are free to lower their total energy by preferentially 
aligning their spins along the magnetic field, leading to a 
temperature-independent Pauli susceptibility. For spin- 
singlet superconductors (both s- and d-wave), the con- 
densate contains an equal number of spin-up and spin- 
down electrons. Therefore, Pauli paramagnetism will 
always favor the normal state over the spin-singlet su- 
perconducting state, and will reduce the superconduct- 
ing critical field H c2 which suppresses superconductiv- 
ity. This effect is called Pauli limiting, with the char- 
acteristic Pauli field Hp determining the upper bound 
of H C 2 3J. Another effect of magnetic field that leads 
to the suppression of superconductivity is orbital lim- 
iting, or suppression of superconductivity when the ki- 
netic energy of the supercurrent around the normal cores 
of the superconducting vortices in Type II superconduc- 
tors becomes greater than the superconducting conden- 
sation energy. The orbital limiting field H® 2 defines H C 2 
in the absence of Pauli limiting. The relative strength of 
Pauli and orbital limiting, the so called Maki parameter 
a = Hfg/Hp, determines the behavior of the system in 
high magnetic field. The prediction of FFLO theory is 
that for a clean Type II superconductor with sufficiently 
large a (for a > 1.8 in the calculations of Ref. 0), a new 
inhomogeneous superconducting FFLO state will appear 
between the normal and the mixed, or vortex, state below 
the critical temperature Tfflo 4]. Within the particu- 
lar realization of Larkin and Ovchinnikov 0, this state 
is characterized by the appearance of a periodic array of 
planes of normal electrons that can take advantage of the 
Pauli susceptibility. 

A number of conventional superconductors were pro- 
posed as candidates for observation of the FFLO state, 



due to their high orbital critical field -ff C 20 and, therefore, 
relatively strong Pauli limiting effect, in the early and 
mid-sixties. Experimental searches, however, yielded null 
results HHSIg. The failure to observe the FFLO state 
was attributed to high spin-orbit scattering rate in these 
compounds jjj In the last decade the FFLO state was 
suggested to exist in heavy fermion UPd2Al 3 (Ref. 
and CeRu 2 (Ref. [HI), based on thermal expansion and 
magnetization data, respectively. Subsequent research 
identified the magnetization feature in CeRu2 as due to 
flux motion 0, and the region of the suggested FFLO 
state in UPd2A L w as shown to be inconsistent with the- 
oretical models 13J . Most notably, multiple phase tran- 
sitions that can be associated with the FFLO state have 
not been observed with a single measurement technique. 

Heavy-fermion superconductor CeCoIns satisfies all re- 
quirements of theory for the formation of the FFLO 
state. It is very clean, with an electronic mean free 
path on the order of microns in the superconducting 
state, which significantly exceeds the superconducting 
correlation length [l^. Its Maki parameter a ~ 3.5 is 
twice the minimum required for the formation of the 
FFLO state ^5|. It was recently discovered that the 
superconducting phase transition changes from second 
to first order at T w 0.3T C for field H || [001], which 
was taken as an indication that Pauli limiting drives the 
physics of CeCoIns at low temperature and high mag- 
netic field 

[HEIl. Th e critical point T was found to be 
in very good agreement with the one predicted by FFLO 
theory for a compound with a = 3.5 Magnetization 
measurements of Tayama et al. [l7| showed that the su- 
perconducting transition in CeCoIns becomes first order 
at a critical temperature To = 0.7 K for both H | [001] 
and H |] [100]. Magnetization measurements of Murphy 
et al. with H || [110] indicated the presence of a sec- 
ond temperature- independent, H w 8 T, anomaly below 
1.4 K, and the authors suggested that these results were 
consistent with FFLO state. 

Materials with quasi- two-dimensional Fermi surfaces, 
which arc likely to exhibit Fermi surface nesting, are ex- 
pected to have more stable FFLO phases when magnetic 
field lies within the 2D-like planes [l9|. De Haas- van 



2 



Alphen studies of CeCoIns revealed that a part of its 
Fermi surface is an undulating cylinder with the axis 
along the (001) direction, characteristic of the quasi- 
two-dimensional systems with planes perpendicular to 
[001] [2(j. These theoretical and experimental ob- 
servations motivated us to perform specific heat investi- 
gation of CeCoIns with magnetic field H _L [001]. 

Specific heat data were collected by employing two 
techniques: the standard quasi-adiabatic method and the 
temperature decay method, where a complete specific 
heat data set for a given field was obtained by differenti- 
ating a single temperature versus time curve, generated 
as the sample was coming into equilibrium with the bath 
starting from high temperature (above 1 K). This tech- 
nique was employed previously to resolve a sharp specific 
heat anomaly associated with the first order supercon- 
ducting phase transition in CeCoIn 5 for H || [001] [l5| . 
and was demonstrated to give high resolution data con- 
sistent with the quasi-adiabatic method. 

Figure ^ shows specific heat data of CeCoIn 5 collected 
with the quasi-adiabatic method, as Sommerfeld coeffi- 
cient 7 = C/T after subtraction of the Schottky anomaly 
tail at low temperature, due to In and Co nuclear lev- 
els Q, for magnetic field H || [110] (panel (a)), and 
specific heat for H || [100] (panel (b)), as a function 
of temperature. The superconducting anomaly at lower 
fields H < 10 T is mean-field-like, with a step in the 
specific heat at T c , similar to the case of H || [001] when 
the field is far from H C 2 [l5| . In this range increasing 
magnetic field simply reduces the magnitude of anomaly, 
without changing the character of the transition. As the 
field is increased further, the trend changes dramatically: 
the magnitude of the anomaly in the specific heat starts 
to increase, and the anomaly itself sharpens up and ac- 
quires symmetric character, characteristic of first order 
phase transitions. The specific heat data indicate that 
the change from second to first order occurs at a criti- 
cal magnetic field Hq ~ 10 T and a critical temperature 
To ~ 1 K. As the superconducting transition temper- 
ature is suppressed by the magnetic field below « 500 
mK, the transition becomes hysteretic (the data for 11.2 
T and 11.4 T in Fig.^b)), proving unambiguously that 
the superconducting transition in CeCoIn 5 at high fields 
close to the critical field H C 2 is indeed first order. At 
a temperature of about 300 mK the specific heat data 
displays an additional anomaly within the superconduct- 
ing state for H > 10 T, which we call a Tfflo anomaly. 
The low temperature region, in the vicinity of the Tfflo 
anomaly, is shown in the insets of Figures QJa) and (b), 
where Tfflo's for different fields are indicated by the 
arrows. The Tfflo anomaly can be described as a step 
followed by a gradual decrease of the specific heat as with 
decreasing temperature, a behavior characteristic of the 
second order phase transition. The Tfflo anomaly is 
observed only in the superconducting state, and disap- 
pears when the superconducting phase transition is sup- 
pressed by magnetic field below TfflOi as illustrated by 
the data for H = 11.4 T in Fig. Ha), or when H < 10 T. 
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FIG. 1: Specific heat of CeCoIn 5 with H 1 [001]. (a) H \\ 
[110] data for fields of 9 T, 10 T, 10.6 T, 11 T, 11.2 T, and 
11.4 T from right to left, collected with heat pulse method. 
Inset: Low temperature region for 10.6 T and 11 T (same 
symbols) emphasizing the Tfflo anomaly. Arrows indicate 
phase transition temperatures from equal area construction, 
(b) H |] [100] Solid symbols: heat pulse data for fields of 9.5 
T, 10 T, 10.5 T, 10.8 T, and 11 T from right to left. Solid 
(dash-dotted) curve is for 11.2 T data collected with the decay 
method with temperature swept up (down) . Dashed (dotted) 
curve is for 11.4 T with temperature swept up (down). Inset: 
Tfflo anomalies for 10.8 T and 11 T. 



The specific heat data collected with the decay method 
for H || [110] is displayed in Fig. [5Ja) as a surface con- 
tour plot in the H — T plane. We can see a clear evo- 
lution of the character of the specific heat anomaly with 
increasing magnetic field from a mean-field-like step to 
a very sharp peak at higher magnetic field, as well as 
the development of the second low temperature Tfflo 
anomaly (a red ridge) in the low temperature/high field 
corner of the H-T plane. By plotting the data as color 
contour plot in Fig. 13b) we can immediately obtain the 
low temprerature/high field part of the phase diagram of 
CeCoIn 5 with H || [110], where both superconducting - 
normal phase boundary T c and the Tp flo anomalies are 
indicated by gray curves. 

The complete H — T phase diagram of CeCoIns based 
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FIG. 2: (a) Electronic specific heat of CeCoIns divided by 
temperature with H || [110] collected with the temperature 
decay method, as a function of field and temperature, (b) 
Contour plot of the data in (a) in the H — T plane. Grey 
lines indicate the superconducting phase transition T c and 
the FFLO-mixed state Tfflo anomaly. Color scale is the 
same in (a) and (b). 



on our specific heat measurements is displayed in Fig. 
for three orientations of the magnetic field, H || [110], 
H || [100] (closed and open symbols in panel (a), re- 
spectively), and H | [001] (panel (b)). The second-to- 
first order change is indicated by To = 1.1 ± 0.1 K for 
H || [110], which is about 10% higher than T for H \\ 
[100]. The To is obtained from the evolution of the spe- 
cific heat anomaly and the magneto-caloric data (not 
shown), with analysis similar to the one performed for 
H || [001] There is anisotropy for the field in the a-b 
plane of CeCoIns 0, . This anisotropy is manifested 
in = 1.1% higher critical field in the [100] direction which 
develops above H = 10 T, the region of the first order 
superconducting transition. The inset (c) of Fig. OJa) 
shows the evolution of the entropy with magnetic field 
H || [100] spanning the region of fields from well into the 
first order (11.4 T) to well into the second order (8.6 T) 
regions of the superconducting phase transition. The en- 
tropy is clearly conserved (all curves collapse on a single 
curve in the normal state (T ss 1.5 K), proving that in 
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FIG. 3: (a): H — T phase diagram of CeCoIn 5 with both H \\ 
[110] (filled symbols) H || [100] (open symbols), (o) and (•) 
indicate the T F FLO anomaly for H || [100] and H || [110], 
respectively. Inset (c): Entropy gain from T — 0.13 K for 
fields of 11.4 T, 11 T, 10.8 T, 10.6 T, 10.22 T, 9.5 T, and 8.6 T 
from left to right. Inset (d): Specific heat jump at the Tfflo 
anomaly obtained from equal area construction, (b) : H — T 
phase diagram for H | [001]. (A) indicate Tfflo anomaly. 
Inset: Sommerfeld coefficient; (o) 5T, solid symbols: 4.9 T, 
4.875 T, 4.85 T, and 4.8 T from left to right. Arrow indicates 
Tfflo anomaly at 4.9 T. Solid lines in (a) and (b) are guides 
to eye for superconducting phase boundaries. 



both regimes the specific heat anomalies are due to the 
same electrons (and no other degrees of freedom) par- 
ticipating in superconducting phase transitions. The in- 
set (d) of Fig. |3fa) shows the magnitude of the step of 
the Tp flo anomaly, obtained via the equal entropy con- 
struction, as a function of magnetic field H | [110]. The 
data are rather linear in field, indicating the tendency of 
the anomaly to disappear for fields less then w 9.9 T. The 
Tfflo anomaly, indicated by solid circles for H \\ [110], 
also appears to extrapolate towards field close to 10 T 
on the ii-axis. The inset in Fig. |3fb) shows low tem- 
perature electronic specific heat (Schottky contribution 
was subtracted) for magnetic field close to H C 2 = 4.95 T 
with H || [001]. The low temperature anomaly Tfflo 
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can also be resolved at 4.9 T, 4.875 T, and 4.85 T. This 
anomaly was not observed for H < 4.8 T. Tfflo = 130 
mK is about half of the value for H || [100]. This in- 
dicates that FFLO state is more stable when magnetic 
field is in the a — b plane of this quasi-2D compound, as 
expected. The tiny high-ficld/low-temperature corner of 
the H — T phase diagram occupied by the FFLO phase 
for if || [001] is indicated by open triangles in Fig.^b). 
The emerging picture therefore is that of a single Tfflo 
phase boundary carving out a high field/low temperature 
part of the superconducting state of CeCoIn 5 . 

A number of theoretical approaches were taken to ex- 
plore the FFLO state, which resulted in a variety of pos- 
sible phase diagrams [H, |22l 0|- O ur data is consis- 
tent with some of these expectations. The first order 
superconducting phase transition for T c < Tq was pre- 
dicted by K. Maki 16] for Type II superconductor with 
strong Pauli limiting. Under these conditions the FFLO 
state was calculated to occur below the same tempera- 
ture To for pure superconductors . Introduction of im- 
purities modifies this picture: the first order normal-to- 
superconducting phase transition is expected to be rather 
insensitive to the impurity scattering, while the FFLO 
state is suppressed to lower temperatures both for the 
s-wave and d-wave [22 pairing. CeCoIns has been 
shown to be a d-wave superconductor in a clean limit |14| , 
with impurity scattering most likely close to the unitary 
limit, based on low temperature thermal conductivity 



measurements. In such Larkin-Ovchinnikov state 

is most likely stabilized in the low temperature/high field 
corner of the superconducting state of the H—T phase di- 
agram , in accord with our data. Recent Monte Carlo 
calculations of the phase diagram of the d x 2_ y 2 supercon- 
ductor in magnetic field 26] indicate that the supercon- 
ducting fluctuations modify the first order phase tran- 
sition below Tq into the nearly discontinuous crossover 
(broadened first order phase transition) , observed exper- 
imentally in CeCoIns. These theoretical considerations 
lead us to conclude that the Tfflo anomaly is indeed 
the vortex state - FFLO state phase boundary. 

In summary, we have observed the low temperature 
specific heat anomaly within the superconducting state 
of CeCoIns in a region of the phase diagram where the 
normal to superconducting phase transition is first order, 
as also demonstrated by the specific heat measurements. 
This transition is conclusively identified as due to the for- 
mation of the spatially inhomogeneous superconducting 
FFLO state, predicted first theoretically about 40 years 
ago. 
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